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ABSTRACT 

We analyse the spatially-resolved stellar populations of 9 local (z < 0.1) Brightest Cluster 
Galaxies (BCGs) observed with VIMOS in IFU mode. Our sample is composed of 7 slow- 
rotating and 2 fast-rotating BCGs. We do not find a connection between stellar kinematics and 
stellar populations in this small sample. The BCGs have shallow metallicity gradients (median 
A[Fe/H] = —0.11 ± 0.1), high central metallicities (median [FdH\[ a / Fe ]=o = 0-13 ± 0.07), 
and a wide range of central ages (from 5 to 15 Gyr). We propose that the reason for this 
is diverse evolutionary paths in BCGs. 67 per cent of the sample (6/9) show ~ 7 Gyr old 
central ages, which reflects an active accretion history, and 33 per cent of the sample (3/9) 
have central ages older than 11 Gyr, which suggest no star formation since z = 2. The BCGs 
show similar central stellar populations and stellar population gradients to early-type galaxies 
of similar mass (M dyn > 10 11,3 M o ) from the ATLAS 3D survey (median [Z/H] = 0.04 ± 
0.07, A[Z/H] = —0.19 d= 0.1). However, massive early-type galaxies from AIhAS 3D have 
consistently old ages (median Age = 12.0 =b 3.8 Gyr). We also analyse the close massive 
companion galaxies of two of the BCGs. These galaxies have similar stellar populations to 
their respective BCGs. 

Key words: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - galaxies: 
evolution - galaxies: kinematics and dynamics - galaxies: stellar content - galaxies: structure 


1 INTRODUCTION 


Brightest Cluster Galaxies (BCGs) are extremely luminous galax¬ 
ies that are usually located in the centre of rich galaxy clusters. 
They have been shown to be distinct from other similarly luminous 
cluster galaxies (e.g. |Hausman & Ostriker] |1978| postman & 
Lauer|[l995l |Lazzati & Chincarini| |1998[ |von der Linden et al. 


2007 


1977 


In the hierarchical scenario of structure formation { [Toomre 


White & Rees 1978 ) galaxies grow in mass and size by 


merging with their neighbours. BCGs are predicted to have a 
more active merger history than lower mass galaxies (White &] 


Rees| 1978, Khochfar & Burkert 20 03} |De Lucia & Blaizot|20071 

Oser et aL||2010[ |Naab et af||2014| ). These galaxies are often 
considered as the extreme end-point of massive galaxy evolution. 


* E-mail:poliva@astro.swin.edu.au 


However, despite being among the most luminous galaxies, and 
generally easy to detect, observations and theory have not reached 
a common point yet, and their evolution is still not fully understood. 


Observations suggest that the mass growth of BCGs evolves 
with time. BCGs accrete their mass at a fast rate until z ~ 0.5, 
thereafter their mass growth slows down ( Lidman et al.|2013]|Lin| 


|et al.|[20l3{ |Qliva- Altamirano et al.||2014[ |Inagaki et al. 2015| >. 

Studies looking at BCG companions have concluded that their 
stellar mass grows through major mergers (^1:3 mass ratios) by 
a factor of 1.8 d= 0.43 at 0.8 < z < 1.5 ( [Burke & Collins|2013| ), 
and mostly by minor mergers (^1:4 mass ratios) by a factor of 
1.1 at z < 0.3 (Ed wards & Patton|2012|>. Major mergers are rare 
at low redshifts, yet still possible (e.g.|Brough et al.|2011||Jimmy| 
|et al.|2013| >. 
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2 P. Oliva-Altamirano et al. 


The recent accretion history of galaxies can be read through 
their stellar population gradients. In the canonical scenario, a 
galaxy’s initial metallicity gradient is set by an initial starburst at 
z ^ 3 and the metallicity decreases in the outskirts, as metallicity 
follows the changes in the gravitational potential ( [Scott et~aT| 
|2009[ |McDermid et al.||2012| ). This gradient can be disrupted by 
violent merging events (major mergers), or reinforced by minor 


mergers ([Kobayashi 

2004, Spolaor et al. 2009, Cooper et al. 2010, 

Pipino et al. 2010). 

Hirschmann et al. (2014) analysed the stellar 


populations of 10 massive halos (10 12 < Mhaio < 10 13 M©) from 
the high-resolution cosmological simulation of |Hirschmann et~aL| 
( |2013| ). They found that major mergers do flatten the metallicity 
gradients. If, as predicted, BCGs have an active merger history, 
including several major mergers, they would be expected to have 
shallower metallicity gradients than lower mass galaxies. However, 
long-slit observations to date suggest that they have a wide range of 
gradients ( [Brough et al.|2007[|Loubser & Sanchez-Blazquez|2012| ). 

Integral Field Unit (IFU) spectroscopy is a valuable tool to 
explore the spatially-resolved kinematics and_stellar populations 
of galaxies. The SAURON (deZeeuw et al. 20021 and ATLAS 313 
(Cappellari et al. 2011) surveys have used IFU spectroscopy to 
explore a significant sample of early-type galaxies in the local 
Universe. [Kuntschner et al.| ( |2010| ) and |McDermid et al.| ( |2015| > 
presented the stellar population analysis of the SAURON and 
ATLAS 3D samples, respectively, finding that 40 per cent of the 
galaxies typically show contributions from young stellar popu¬ 
lations connected to low mass, fast rotator systems. In contrast, 
they find that slow rotators are generally consistent with old (^10 
Gyr) stellar populations. The most massive systems (stellar mass 
^ 10 10 5 M 0 ) have the flattest metallicity gradients. However, the 
ATLAS 3D sample contains 21 galaxies with dynamical masses 
greater than 1O 113 M 0 , and only one of those is a BCG (M87). 


| Jimmy et al.| ( |2013| ) analysed the kinematics and photometry 
of a sample of 10 BCGs observed with the VIMOS IFU, of which 
4 have close massive companions. If BCGs were the product of 
many minor mergers they would be expected to be slow-rotating 
galaxies. [Jimmy et al.| ( |2013| found that 30 per cent of the BCGs 
in their sample are fast rotators. The simulations of |Naab et al.| 
< |2014| ), predict that angular momentum mostly depends on the gas 
content of the galaxies involved in the interaction, suggesting that 
the slow or fast rotation could be a temporary state. |Jimmy et al. 
(|2013| ), also find through photometric analysis (G — M 20 ; Lot z 


|et al.| 2008) that 40 per cent of the galaxies in the sample have 
undergone a minor merger within the last 0.2 Gyr. 


Our paper is the third in a series analysing the spatially-resolved 
spectra of BCGs (after |Brough et al.||201 1~[ | Jimmy et al.||2013| ). 
We present here the spatially—resolved stellar populations of the 
BCG sample presented in | Jimmy et al.| |2013| ). This is the first 
IFU analysis dedicated to the stellar populations of BCGs. We will 
investigate whether BCG accretion histories are different from 
those of the general massive early-type galaxy population by com¬ 
paring our measurements to those from SAURON and ALTLAS 315 . 


We present the sample selection and observations in Section 
2. Section 3 describes the stellar kinematic and photometric 
measurements from | Jimmy et al.| ( |2013|. It also describ es the 
early-type galaxy samples of SAURON (Kuntsch ner et al.|[2010| > 
and ATLAS 3D (Mc Dermid et al.|2015 1 to which we will compare 
our sample BCGs throughout the paper. Section 4 presents the 


stellar population analysis. Section 5 summarises the main results. 
These are later discussed in Section 6. Our conclusions are 
presented in Section 7. The cosmology adopted throughout this 
paper is Ho=70 kms -1 Mpc -1 , Om=0.3, Q a =0.7. 


2 OBSERVATIONS 

In this section we summarise the observations made and the data re¬ 
duction process. These are described in full in |Jimmy et al.|f20T3) . 


2.1 Spectroscopic Measurements 

The BCG sample is selected from [von der Linden et al.| ( |2007| ). 
The galaxies are part of the C4 cluster catalogue (Miller 2005) 
of the third data release of the Sloan Digital Sky Survey (SDSS; 

| York et ak||2000| . These observations consist of 10 BCGs, 4 of 
which have massive companions within ~ 10 ,/ (corresponding to 
~ 18 kpc at z < 0.1). We use the same nomenclature as Jimmy] 
|et~aL| (2013), i.e. we present each cluster as the last 4 digits in 
the SDSS flag, rather than SDSS-C4-DR3 number. The galax¬ 
ies were observed with the Very Large Telescope using the IFU 
mode of the VIMOS spectrograph ( |Le Fevre et al.|[2003] >, with 
the high-resolution blue grism, which has a spectral resolution of 
0.51 A/pixel. The observations were made in two sets, April to Au¬ 
gust of 2008 and April to July of 2011 (Prog. ID 381.B-0728 and 
Prog. ID 087.B-0366, respectively). The galaxies used in this study 
have a spatial sampling of O.GT'Vpixel with a field-of-view (FOV) 
of 27" x 27". The average seeing of the observations is 0.9". The 
rest wavelength range is ~ 3900 to 5600 A. 


2.2 IFU Data Reduction 

The IFU data reduction consists of two different stages: (a) The VI¬ 
MOS Pipeline ( |Izzo et al.|2004) which generates the calibrations 
files (fibre identification, master bias, etc.), and does a first order 
flux calibration. This flux calibration corrects the spectrum shape, 
using the standard stars observed on the same night as the galax¬ 
ies. The VIMOS FOV is formed by four quadrants. The VIMOS 
pipeline reduces each quadrant separately. As a result we obtain 
the science spectrum for each spaxel in each quadrant, (b) We use 
our own IDL routines to mask the bad fibres, subtract the sky for 
each quadrant and then combine the quadrants into a three dimen¬ 
sional data cube. The multiple exposures for each observation are 
combined using a 5 a clipped median. This code is publicly avail¬ 
ably Finally, we flux calibrate the spectrum using the photometric 
standard stars. 


3 PREVIOUS MEASUREMENTS 

| Jimmy et al.| ( |2013] ) measured the stellar kinematics and photome¬ 
try of the galaxies in our sample, we summarise their method and 
results in this section. 


1 http://galaxies.physics.tamu.edU/index.php/Jimmy#Code 
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Stellar Populations of Brightest Cluster Galaxies 3 


3.1 Stellar Kinematics 

First a signal-to-noise ratio (SNR) cut of 5 across all the spaxels 
(spatial pixels) was applied. The spaxels were then re-binned to a 
minimum SNR of 10, using the spatial binning Voronoi code of 
Cappellari & Copin (2003). The velocity and line-of-sight velocity 
dispersion were computed using the penalised fitting scheme of 
|Cappellari & Emsellem| (pPXF; |2004| ), and the MILES (Medium- 
resolution Isaac Newton Telescope Library of Empirical Spectra; 
|Sanchez-Blazquez et al.||2006| ) library stellar templates. pPXF 
fits the stellar library templates to the absorption line features of 
the BCG spectra, giving the redshifts and the broadening of the 
spectral lines. 


The angular momentum was characterised by the Ar pa¬ 
rameter defined by Emselle m et al.| ([2007). It is calculated as 
follows: 


Ar 


(R\v\) 

(R^V 


( 1 ) 


where R represents the radius of the galaxy, V is the stellar ve¬ 
locity and a the velocity dispersion. The numerator and denom¬ 
inator are luminosity weighted. A higher Ar represents a higher 
angular momentum. The ellipticity (e) at the effective radius of 
each galaxy was measured using the publicly available IDL routine 
find_galaxy .pro developed by Michele CappellarQ Following 
|Emsellem et al.| ( [T011| ), the values of Ar and e can be used to dis¬ 
tinguish fast and slow rotators (FR and SR respectively) by using 
the threshold: 


Ar ^ (0.31 =b 0.01) x v/i, (2) 

where FRs lie above this threshold and SRs lie below. 


The dynamical mass was measured using the standard equation 
given in Cappellari et al. ( 2006). 


Mdyn — 


5 R e al 


(3) 


where a e is the aperture corrected velocity dispersion of the inte¬ 
grated spectrum within the effective radius; G is the gravitational 
constant. In Table[l]we summarise the relevant kinematic results. 7 
of the BCGs are SR and 2 are FR. For cr e , e e and A R e we refer the 
reader to Table 2 of |Jimmy et al.H2013) . 


3.2 Photometric Analysis 

| Jimmy et al.[ ( |2013] ) analysed the photometry of this sample using 
images from SDSS Data Release 3. They measured the effective ra¬ 
dius by fitting a 2D de Vaucouleurs profile. They also analysed the 
presence of recent mergers using the Gini, and M 20 coefficients 
( |Lotz et aL||2008| >. This method studies the distribution of light 
looking for irregularities that could indicate morphological signa¬ 
tures of mergers. A galaxy is a merger candidate if it crosses the 
threshold: 


G ^ —0.14M2O + 0.33 (4) 

Where M 20 is the 2nd order moment of the brightest 20 per cent of 
pixels, and G is the Gini coefficient. In the case of gas-poor galax¬ 
ies like our sample, a galaxy will be above the threshold if it is 
currently merging or has merged in the last 0.2 Gyr ( |Lotz et ak| 


Table 1 . Kinematic properties of BCGs and their companions 
from |Jimmy et al.| < j2013| . Seven of the BCGs are slow rotating 
(SR) and two are fast rotating (FR). Four of the BCGs show 
photometric signs of merging. 


Galaxy 

z 

log M dyn M© 
± 0.01 

M") 

±0.01 

Merging? 

G - m 20 

FR/SR 

BCGs 

1027A 

0.090 

11.79 

6.98 

y 

SR 

1042 

0.094 

11.83 

7.22 

n 

SR 

1050 

0.072 

11.78 

8.43 

n 

SR 

1066 

0.083 

11.62 

5.07 

y 

SR 

2001 

0.041 

11.38 

5.84 

n 

SR 

2039 

0.082 

11.86 

8.82 

n 

SR 

2086 

0.083 

11.60 

4.83 

y 

SR 

1048A 

0.077 

11.59 

5.17 

y 

FR 

1261 

0.037 

11.32 

5.76 

n 

FR 

Comp 

1027B 

0.090 

11.17 

4.39 

y 

FR 

1048B 

0.080 

10.51 

1.08 

y 

FR 

1048C 

0.074 

10.54 

1.24 

y 

FR 


I 20 TT) . 4 of the BCGs in the sample are merging. The relevant kine¬ 
matic and photometric results of | Jimmy et al.| ( [20T3] ) are presented 
in Table [T] 


3.3 Stellar populations from the SAURON and ATLAS 315 
samples 

Throughout the paper we compare our observations to those of 
early-type galaxies of similar mass observed by ATLAS 315 (which 
includes the SAURON galaxy sample). The ATLAS 315 sample 
is composed of 260 field and cluster early-type galaxies and 
it only contains one BCG, M87. The spatially-resolved stellar 
populations (central values and gradients) of the SAURON sample 
were presented in Kuntschne r et al.| ( |2010| . The ce ntral stellar 


populations of the ATLAS 315 sample were presented in 
|et al.|p0l5) . We therefore compare our central stellar populations 
with the whole ATLAS 315 sample and the stellar population 
gradients only with the SAURON sample. 


McDermid 


Kuntschn er et al.| ( |2010| > and |McDermid et al.| < |2015| ) use the 
stellar models of |Schiavon| ( |2007| > in the Lick/IDS system 
jWorthey & Ottaviani 1997) to measure the stellar population 
parameters of age, total metallicity [Z/H] and abundance of alpha 
elements [a7Fe]. For our spectral fitting, we use models with 
only solar abundances, i.e. [a7Fe] = 0. In this case [Fe/H] is 
effectively a measure of the total metallicit)[^] Therefore, we 
directly compare the ATLAS 315 total metallicities to our measured 
metallicities throughout the paper. The central stellar populations 
in the ATLAS 315 sample correspond to an aperture of 0.125 R e . 


We compare the median, and standard deviation for the two 
samples in the same mass range (M d yn > 10 n ' 3 M©). In order to 
compare our BCGs with a non-BCG early-type galaxy sample, we 
do not include the BCG M87. This gives a comparison sample of 


2 http://www-astro.physics.ox.ac.uk/~mxc/idl/ 


[Fe/H] = [Z/H] - 0.75*[a7Fe] j Conroy & van Dokkum 


2012 1 
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20 massive early-type galaxies from the ATLAS 30 sample. We 
highlight M87 as a blue-filled circle in the figures. 


4 STELLAR POPULATION ANALYSIS 

The study of stellar populations requires high SNR spectra. In 
order to secure a high enough SNR we use our own python routing 
to stack the spaxels within annuli for each galaxy. To identify 
the annuli we follow the total flux in the wavelength—integrated 
galaxy image. This indirectly maps the galaxy morphology, i.e. 
the shape of the annuli is determined by the galaxy’s morphology. 
Each spectrum is shifted to rest frame wavelength before stacking 
using the velocity measurements obtained from pPXF. The spectra 
are then broadened to a reference velocity dispersion, a (the 
maximum velocity dispersion of the galaxy). This reduces the 
dilution of the spectral features due to rotational broadening and 
allows us to have a fixed and constant velocity dispersion when 
measuring the stellar population parameters. As a result, we have 
one spectrum per annulus per galaxy. 

From the 10 BCGs and 4 companion galaxies presented in 
|Jimmy et al. 1(201 3) , the BCG 1153, and the companion of BCG 
1066 have too low SNR for stellar population analysis. Our final 
sample thereby consists of 9 BCGs, 4 of them with close massive 
companions. For 2 of these (1027, 1048) it was possible to resolve 
the companions as well. In those cases we refer to the main galaxy 
as 1027A and 1048A, and the companions as 1027B, 1048B, and 
1048C. 

In Fig [T] we show the annular distribution of 3 representative 
galaxies. The upper panels are the flux-collapsed YIMOS image. 
The lower panels are the annular distribution per galaxy. The 
annuli cover up to 1 R e in each galaxy. The central aperture has 
been defined as 0.2 d= 0.03 R e , allowing the central annulus to 
contain two or more spaxels (for most of the galaxies). The median 
value of the SNR of the final sample stacked spectra is ~ 35 A -1 
with the majority having SNR > 20 A -1 . Only the outermost 
annuli in the companion galaxies 1048B and 1048C drop below 
this threshold (SNR = 15 A -1 ). 


4.1 Stellar Population Models 


We use the stellar population models of Vazdekis et al. 


based on the MILES stellar library ( [Sanchez-Blazquez et a . 


( 20101 , 


2006), 


to estimate the ages and metallicities of galaxies in our sample. 
The|Yazdekis et al.] ( |2010| ) models use the Padova 2000 ( [Girardi 
et al.|2000| ) isochrones which cover a metallicity range of [Fe/H] 


= [-1.7, 0.4] and ages [2.1 x 10 7 ,1.7 x 10 10 ] yr. After testing our 
spectra with the whole age range, we restrict our analysis to single 
stellar population sequences of [1 x 10 9 ,1.7 x 10 10 ] yr, divided 
into 40 age bins. 


We also considered the stellar population models of PEGASE-HR 
(|Le Borgne et al.||2004J coupled to the ELODIE (|Prugniel et al.| 
[2007 ) stellar library. These models are known for their high 
spectral resolution (0.55 A) and for allowing the user to choose 
the initial mass function along with other physical ingredients. Its 
wavelength coverage is consistent with our data, (A > 3900 A). 
Furthermore, the stellar library has the same age range as MILES. 
We find both libraries to give similar results. Both metallicities 
and ages are consistent within lcr uncertainties. However, the 
flux calibration of the MILES library over a wide spectral range 
is more appropriate for our spectral fitting of unresolved stellar 
populations. The average spectral resolution of MILES (2.3 A) 
is also closer to that of the VIMOS observations (2.1 A). This 
minimises information loss when broadening the library spectra. 
We therefore use the [Vazdekis et al.| ( |2010| ) models in the results 
presented here. 


4.2 Full Spectrum Fitting 

To estimate the stellar population parameters we used the full 
spectrum fitting technique (e.g. |Koleva et al.|2008| ). This technique 
features some advantages over classical methods, i.e. the Lick/IDS 
system (e.g. |Worthey & Ottaviani||1997| ), as it exploits all the 
information contained in the spectra, pixel by pixel, independently 
of the spectrum shape. It also allows analysis at medium spectral 
resolutions (< 3 A), in contrast to the Lick/IDS system which has 
a low spectral resolution, (> 8 A). We use the STEllar Content and 
Kinematics via Maximum A Posteriori algorithm (STECKMAP; 
|Ocvirk et al.|2006afb} to extract the ages and metallicities from our 
spectra. 


The analysed rest wavelength range (3900-5600 A) covers 
the absorption line indices: Can (Ca II 3933), CaH (Ca II 3968), 
H 6, H 7 , Ca G (Ca I 4307), H£, Fe5015, Mgb5175, Fe5270. Our 
age and metallicity measurements are based on the full spectrum. 
However, we mask the weaker B aimer indices H 5, and H 7 to 
ensure a clean comparison with the ATLAS 3D data (which age 
measurements are based on H /3). The enhanced lines Ca k, CaH, 
and Ca G provide a strong case for the existence of old stellar pop¬ 
ulations. Furthermore, where there are mixed stellar populations, 
the ratio between these lines indicates the relative importance of 
the young stellar populations ( [Sanchez Almeida et al.[2012] ). 

In the following subsections we describe the method we im¬ 
plement to estimate the metallicities and ages. 


STECKMAP uses Bayesian statistics to estimate the stellar 
content of the spectra. It is based on a non-parametric formalism. 
The code returns a luminosity-weighted age and metallicity 
distribution. This comes from flux-normalising the stellar library 
(rather than mass normalisation). The results are a proxy to the 
star formation history of the galaxy. The method is regularised 
by a Laplacian kernel in order to avoid chaotic oscillations. To 
prevent systematic errors from poor flux calibration, the code 
produces a non-parametric transmission curve which represents 
the instrumental response multiplied by the interstellar extinction. 
We further mask the spectral region around the weak emission 
lines so that they do not interfere with the stellar population model 
fitting (e.g, [Neill] 3868.71, [OIII] 4959, 5007, [NI] 5198, 5200). 

In order to obtain accurate, robust results we first test the 
spectrum by measuring the stellar population parameters using two 


4 http://astronomy.swin.edu.au/~poliva/codes/annuli_stacking/spectra.py 


5 [Fe/H] = logio (Z/Z©), and Z© 


0.0189 


( Ande rs & Grevessejl 989i 
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Figure 1. Three representative BCGs: 1050 (SR), 2039 (SR), 1261 (FR). Upper panels: The flux-collapsed VIMOS image (27" x 27"). Lower panels: The 
annular distribution for each galaxy. The central annulus, corresponds to the aperture 0.2 d= 03 R e . The annuli all extend to 1 R e . 


different age initial conditions: (a) a flat stellar age distribution, 
(b) a random Gaussian distribution of ages (see |Ocvirk||2QTT) . 
We then compare the results from both runs expecting them to be 
consistent. If this test is successful, we proceed to calculate the 
final value of the stellar population parameters. The final estimated 
luminosity-weighted age and metallicity are the median values 
from 150 Monte Carlo realisations. In each Monte Carlo realisation 
the initial condition is a random Gaussian age distribution that 
is later refined through iterations until it reaches the best fit. The 
measurement uncertainties are the standard deviations of the 150 
Monte Carlo realisations. In Fig [2] we show the spectrum of the 
central annulus of 2039. The black line is our data, the red line 
represents the best fit to the data, the vertical dotted lines indicate 
the significant lines. 

Unfortunately, STECKMAP, in conjunction with the 
MILES/ELODIE stellar libraries, has the disadvantage of be¬ 
ing tied to solar abundance ratios ([a/Fe] = 0). This means that 
were not include very high metallicities and a-element [a/Fe] 
abundance in this model. BCGs are known to have high metallici¬ 
ties and super solar [a/Fe] ratios (e.g. |von der Linden et al.|2007) . 
We therefore explore the impact of this on our fits using a Lick 
index analysis. 

4.3 Lick/IDS System 

The stellar population parameters can also be estimated using 
the Lick/IDS system (e.g. [Worthey & Ottaviani| |1997| . This 
method measures the equivalent width of the absorption features 
at a fixed IDS resolution (> 9 A) to later compare them with 
stellar population models that provide the corresponding age and 
metal licity values (e.g.|Worthey & Ottaviani|1997l |Vazdekis et al.| 
1 1997 1 proctor & Sansom|2002||Schiavon|2007|>. 



4000 4200 4400 4600 4800 5000 5200 5400 
wavelength [A] 


Figure 2. Rest wavelength spectrum of the central annulus of BCG 2039. 
The black line represents our data, the red line shows the best fitting stellar 
populations from STECKMAP. The green line shows the residuals of the 
fit. The shaded regions show the regions that were masked on the fit. The 
major features are marked with dotted vertical lines. The recovered stellar 
populations are [Fe/H][ a / Fe ] =0 = 0.13 =b 0.01 and Age = 9.6 =b 1.0 Gyr. 
The fit is a good match to the data. 


To probe the robustness of our results obtained using STECKMAP 
we test our VIMOS spectra with the Lick/IDS system (see also 
Appendix A). We use the same stellar models and absorption lines 
as the ATLAS 315 team to measure the [a/Fe] abundances, ages 
and metallicities i.e. |Schiavon| (2007) , H/3, Fe5010, and Mgb. 
In Fig [3] we show the central [a/Fe] of our BCGs (green-filled 
and red-open squares) and the ATLAS 315 sample (blue crosses) 
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log M dyn (M@) 

Figure 3. Central a-abundance [a/Fe] as a function of galaxy dynam¬ 
ical mass. The blue crosses r epresent the early-type galaxies from the 
ATLAS 30 jMcDermid et al. ||2015| survey with M87 shown as a blue- 
filled circle. The slow rotator BCGs are shown as green-filled squares. 
The fast rotator BCGs are shown as red-open squares. The dashed line 
shows the mass range (M^ n > lO 11 -^©) used in the comparison be¬ 
tween the two samples. We find that the BCGs have similar a-enhanced 
ratios (median [a/Fe] = 0.17 ± 0.04) to the early-type galaxies (median 
[a/Fe] = 0.24 ± 0.03). 

as a function of mass. We find the BCG [a/Fe] values (median 
[a/Fe] = 0.17 =b 0.04) to be consistent with those of the ATLAS 315 
massive galaxies (median [a/Fe] = 0.24 zb 0.03). 

In Fig [4] we show the estimated metallicities and ages from 
both methods, STECKMAP and Lick indices. The Lick 
ages (median Age = 10.0 zb 1.1 Gyr) and metallicities 

(median [Z/H] = 0.18 zb 0.02) are consistent with those 
from STECKMAP (median Age = 8.9 zb 3.3 Gyr, median 

[Fe/H][ a / Fe ]=o = 0.13 zb 0.07) within la error. Our STECKMAP 
results are also consistent with the stellar populations analysed by 
|Gallazzi et al.| ( |2005| ) from the SDSS DR4 spectra (of the galaxies 
in our sample). 

We adopt as the final result the stellar population parameters 
measured with STECKMAP. We discuss these results in the 
following sections. 


5 RESULTS 

We present here the main results from our BCG stellar population 
analysis. 

5.1 Central Stellar Populations 

Our stellar population measurements come from a luminosity- 
weighed distribution, therefore our results are sensitive to the 
brightest and youngest population of stars in the galaxy (e.g. 
jTrager & Somerville||2009| . The central values are measured in 
an aperture of 0.2 zb 0.03 R e and are given in Table [2] We find 
that 6 out of 9 BCGs (67 per cent of the sample) have central 
intermediate ages (5 Gyr < Age <10 Gyr), and 3 out of 9 (33 
per cent of the sample) are centrally old (Age >11 Gyr). The 
median central age is 8.9 zb 3.3 Gyr. The BCGs have homogeneous 



[Z/H ] Uck 



Figure 4. Comparison between STECKMAP and Lick estimated stellar 
populations. The green-filled squares represent the central stellar popula¬ 
tions of the SR BCGs and the red-open squares represent the FR BCGs. 
The line shows the one-to-one relationship. The [Fe/H] [ a /Fe]=o is equiv¬ 
alent to [Z/H] (see section 3.3). Our STECKMAP and Lick measurements 
agree within la error. 

super-solar metallicities (median [Fe/H][ a / Fe ] =0 = 0.13 zb 0.07) 
in their central regions. 

Fig [5] shows the central metallicities (upper panel) and ages 
(lower panel) as a function of galaxy mass. The SR and FR BCGs 
are presented as green-filled squares and red-open squares respec¬ 
tively. The companion galaxies (all FRs) are shown as open stars. 
We do not find any significant difference between the SRs and FRs, 
consistent with McDeri md et al.| 42015]) . The crosses represent 
the central stellar populations of the ear ly-type galaxies in the 
ATLAS 315 sample (McDermid et al. 2015). From a Kolmogorov - 
Smirnov (K-S) test we find that the hypothesis that the massive 
early-type galaxy metallicities and BCG metallicities come 
from the same distribution, can not be rejected at a 10 per cent 
confidence level (P-value 0.11). The BCGs have similar central 
metallicities (median [Fe/H][ a / Fe ] =0 = 0.13 zb 0.07) to early-type 
galaxies in ATLAS 315 (median [Z/H] = 0.04 zb 0.07) at fixed mass. 

The BCGs have slightly younger central ages (median 
Age = 8.9 zb 3.3 Gyr) compared to the ATLAS 315 galaxies 
at the same mass (median Age = 12.0 zb 3.8 Gyr). From a K-S test 
we find that the age distribution of early-type galaxies and BCGs 
are different (P-value 1.2 x 10 -6 ). We also find that the BCGs 
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Figure 5. Central stellar populations as a function of dynamical galaxy mass. The SR BCGs are shown as green-filled squares. The FR BCGs are shown as 
red-open squares. The FR companion galaxies are shown as open stars. The blue crosses represent the central stellar populations of the early-type galaxies in 
the ATLAS 315 sample. The blue-filled circle is the BCG M87 from the ATLAS 315 sample. The dashed line shows the mass range (M d yn > 10 11-3 M©) used 
in the comparison between the two samples. Upper panel: Central metallicities as a function of dynamical galaxy mass. The BCGs show similar metallicities 
(median [Fe/H] [ OL /Fe}= o = 0.13 d= 0.07) compared to the ATLAS 315 galaxies at the same mass (median [Z/H ]= 0.04 d= 0.07). Lower panel: Central ages 
as a function of mass. The BCGs central ages (median Age = 8.9 =L 3.3) are consistent with the ages of the ATLAS 315 galaxies of the same mass (median 
Age = 12.2 zb 1.3) within 2cr error. The central stellar populations of BCGs show little scatter or dependence on their angular momentum. 


Table 2. Central stellar populations of BCGs and their compan¬ 
ions. 



Galaxy 

log [Fe/H] 

log [Fe/H] 

Age 

Age error 




error 

(Gyr) 

(Gyr) 

BCG SR 

1027A 

0.18 

0.049 

8.9 

LI 


1042 

0.17 

0.049 

6.4 

1.1 


1050 

-0.003 

0.032 

16.4 

1.1 


1066 

0.235 

0.021 

5.8 

1.1 


2001 

0.163 

0.097 

12.6 

1.3 


2039 

0.132 

0.039 

9.6 

U 


2086 

0.02 

0.068 

13.2 

1.2 

BCG FR 

1048A 

0.121 

0.069 

7.1 

1.2 


1261 

0.084 

0.029 

8.8 

LI 

BCG 

median 

0.132 

0.048 

8.9 

1.1 

Comp 

1027B 

0.127 

0.08 

5.4 

1.2 


1048B 

-0.203 

0.139 

6.6 

1.3 


1048C 

0.198 

0.069 

7.6 

1.2 

Comp 

median 

0.126 

0.079 

6.5 

1.2 


NOTE: The errors listed are the measurement uncertainties. 


have similar central stellar populations to the companion galaxies. 

In summary, our sample BCGs show homogeneous central 
metallicities and a wide range of ages. They show similar 
metallicities to other early-type galaxies of similar mass. 


Table 3. Stellar population gradients of BCGs and their com¬ 
panions. 



Galaxy 

A [Fe/H] 

A [Fe/H] 

A Age 

A Age 




error 


error 

BCG SR 

1027A 

-0.115 

0.08 

0.004 

0.061 


1042 

0.028 

0.094 

0.019 

0.162 


1050 

0.092 

0.042 

-0.197 

0.066 


1066 

-0.414 

0.158 

0.054 

0.116 


2001 

-0.054 

0.225 

-0.134 

0.092 


2039 

-0.171 

0.09 

0.03 

0.036 


2086 

0.05 

0.026 

-0.091 

0.102 

BCG FR 

1048A 

-0.213 

0.221 

0.135 

0.119 


1261 

-0.198 

0.096 

0.034 

0.045 

BCG 

median 

-0.115 

0.094 

0.019 

0.092 

Comp 

1027B 

-0.17 

0.006 

0.214 

0.089 


1048B 

0.174 

0.35 

0.066 

0.061 


1048C 

-0.11 

0.163 

0.133 

0.125 

Comp 

median 

-0.110 

0.160 

0.133 

0.088 


NOTE: The errors listed are the measurement uncertainties. 

5.2 Age and Metallicity Profiles 

The stellar population gradients and their uncertainties are 
measured using a linear log-log chi-squared fitting routine. The 
uncertainties are the standard deviation of the fit. We summarise 
the stellar population gradients (A [Fe/H] and A Age) in Table [3] 
The metallicity and age profiles for the BCGs and companion 
galaxies are illustrated in Appendix B. Most of the BCGs in our 
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Figure 6. Stellar population profiles of BCGs and early-type galaxies of 
similar mass. The solid lines represent the median stellar population gra¬ 
dients of the BCGs in our sample. The dashed lines represent the median 
stellar population gradients of the most massive galaxies in the SAURON 
sample. The grey (BCGs) and blue (SAURON) shaded regions indicate the 
error on the median. These values are specified in each panel. BCGs and 
early-type galaxies at the same mass have similar gradients. 


sample have shallow metallicity gradients, A[Fe/H] > —0.3, 
except for 1066 which has a gradient of A[Fe/H] —0.41 d= 0.1. 
The median value is A[Fe/H] = —0.11 ±0.1. The companion 
galaxies also show shallow metallicity gradients. However 1048B 
and 1048C have R e values close to the seeing FWHM, which may 
act to dilute the measured gradient in these galaxies. 

We summarise the stellar population gradients we observe 
for the BCGs in Fig [6] and compare that to the gradients of the 
SAURON galaxies at the same mass range. The offset between the 
two profiles illustrates the differences between the central stellar 
populations in the two samples. 

In Fig [7] we show the metallicity (upper panel) and age (lower 
panel) gradients as a function of galaxy mass. The BCGs have shal¬ 
low metallicity (median A[Fe/H] = — 0.11±0.1) and age gradients 
(median AAge = 0.02 ± 0.03), similar to those of the SAURON 
early-type galaxies at the same mass (A[Fe/H] = —0.19 ±0.1, 
AAge = 0.04 ± 0.05). We do not find any correlation between the 
stellar kinematics and the stellar population gradients. 


6 DISCUSSION 

BCGs represent the extremely massive end of the early-type galaxy 
population. These galaxies live in high-density environments com¬ 
monly surrounded by many companions. We have presented here 
the first integral field analysis of the radial stellar populations of 9 
BCGs up to 1 R e . 


6.1 Stellar Ages 

Hydrodynamical simulations of early-type galaxies (in less dense 
environments than BCGs) predict more massive galaxies to be 
older than less massive galaxies, such that at masses > 10 10 ' 5 M 0 
the galaxies are older than 10 Gyr (e.g. Naab et al. 2014, Peeples 


|et al.||2014| jHirschmann et al.||2013| ) and show passive evolution 
from z = 2. We find that 3 out of 9 BCGs in our sample have old 
central ages (> 12 Gyr), in agreement with this prediction. These 
3 galaxies are also consistent with the massive early-type galaxies 
from ATLAS 3 ^ (median Age — 12.0 ± 3.8 Gyr). 


However, 6 out 9 BCGs in our sample have central interme¬ 
diate ages (5 Gyr ^ Age <10 Gyr). Previous observations have 
shown that these intermediate ages in BCGs are not unusual. 
|Loubser et al.| ( |2009[ ) analysed a large sample of 49 BCGs, and 
found that 24 of them (49 per cent of the sample) are younger than 
9 Gyr old. Fitzpatrick & Graves ( 2014) found that, at fixed velocity 
dispersion and surface brightness, central galaxies in SDSS are 
younger than satellite galaxies. La Barbera et al. (2014) also found 
that central galaxies have younger ages and higher metallicities 
than isolated early-type galaxies. Consistently, these 6 BCGs sit 
at the younger age limit of the massive early-type galaxies from 
ATLAS 3 ^ (Fig[5J. 


Many hypotheses have tried to explain why such massive 
galaxies as BCGs have intermediate age stellar populations. A 
common prediction is that gas cooling from the intra-cluster 
medium may be forming stars which will result in young ages 
(e.g. |Edwards et al.pOOTllQ’Dea et al.|2~008l|Bildfell et al.|2008[ 


|Loubser|2014| >. For most of the clusters considered here, X-ray 

imaging is not available, so that we cannot assess whether they host 
a cool core or not. Only two of our galaxies have been confirmed 
to be hosted by a non-cooling flow cluster ( [White et al.|1997| and 
for those we find contradictory results: 2001 is one of the oldest 
galaxies in our sample, consistent with the cool core hypotheses, 
however, 1066 shows a young central region and a positive age 
gradient, suggesting that cool cores are not the explanation for the 
intermediate ages observed. 


The semi-analytical model of Tonini et al.] ( |2012| ) focusses 
on BCGs rather than on the general population of early-type 
galaxies, predicting that BCGs have more prolonged star formation 
as a result of their active merger history. The BCGs experience 
continuous bursts of star formation across cosmic time, generating 
many stellar populations super imposed on one another. Our 
age measurements are luminosity-weighted, which means they 
reflect the youngest stellar population of the galaxy. Therefore, the 
intermediate ages we find suggest that the last star formation event 
could have taken place at z ~ 1 when galaxy mergers are more 
likely to be gas rich. 


6.2 Metallicities 


The BCGs in our sample show very homogenous central metal¬ 
licities (median [Fe/H] [ a /Fe]=o — 0-13 ± 0.07). These high 
metallicities are consistent with previous long-slit and fibre 
observations of BCGs (Broug h et al.|2007l von der Linden et al.| 


|2007| [Loubser et al.|2009||Eigenthaler & Zeilinger|2013| ) and are 

in agreement with the hypothesis of continuous star formation 
events at high redshifts. 


We find that the BCGs we observe have a range of metallicity 
gradients, from flat to shallow (median A [Fe/H] — —0.11 ± 0.1) 
similar to other massive early-type galaxies at similar mass 
(A[Z/H] = -0.19 ± 0.1; Fig[7). 

Hydrodynamical simulations predict that galaxies that form 
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Figure 7. Metallicity (upper panel) and age (lower panel) gradients as a function of dynamical galaxy mass. The SR BCGs are shown as green-solid squares. 
The FR BCGs are shown as red-open squares. The blue dots represent the galaxies in the SAURON sample. The blue-filled circle is the BCG M87 from the 
SAURON sample. The dotted line indicates a flat gradient. The dashed line shows the mass range (Md yn > 10 1:L * 3 Mq) used in the comparison between 
the two samples. The BCGs (A[Fe/H] = —0.11 db 0.1) and the most massive early-type galaxies (A[Fe/H] = —0.19 db 0.1) have similar stellar population 
gradients. 


through dissipative core collapse have typical metallicity gradients 
A[Fe/H] ~ —0.4 (e.g |Kobayashi|2004[ [Hirschmann et al.|2014| . 
This is significantly steeper than the gradients we observe. How¬ 
ever, simulations also show that this initial gradient can later be 
affected by accretion of external stellar populations, i.e. mergers 
( |Hirschmann et al.|20131|Martizzi et al.|2014| ). 

Dissipationless major mergers make the stars lose their or¬ 
bits and move randomly within the distribution of the galaxy, 
inducing stellar population mixing and flattening the gradients 
( |Hopkins et al.|[2009] >. This suggest that the BCGs as well as 
the massive SAURON galaxies have gone through at least one 
recent dissipationless major merger since z < 1 fKobayashi|2004[ 
|Hirschmann et al.|2014| ). 

6.3 Merger Histories 

From our analysis we conclude that BCGs have diverse evolution¬ 
ary paths. 3 out of 9 BCGs in our sample show old and metal-rich 
central stellar populations, and shallow metallicity gradients. This 
suggest that their stars were formed in-situ at z > 2. Thereafter 
the galaxies grow in mass and size by at least one major merger 
and many minor mergers (e.g. Kobayashi 2004, Hirschmann 
|erah1|2014| . These galaxies are similar to the massive early-type 
galaxies in the SAURON and ATLAS 3i:> sample which also have 
old central stellar populations and shallow metallicity gradients. 

The rest of the sample (6 out of 9) BCGs have intermediate 
central ages, high central metallicities, and shallow metallicity 
gradients in BCGs. This implies that these galaxies have ex¬ 
perienced active accretion histories throughout cosmic time, as 
predicted by semi-analytical and dark matter simulations |De] 
|Lucia & Blaizot|2007||Tonini et al.|2012||Laporte et al.|2013] >. The 

dense environment where BCGs evolve allows them to experience 


many mergers. These mergers will trigger star formation at high 
redshifts, and will disrupt the metallicity gradients at z < 1, given 
that the fraction of gas in the merging galaxies decreases with time. 

| Jimmy et al.| ( [2013) found that 4 of the 9 BCGs studied here 
show photometric signatures of minor mergers. The effect of minor 
mergers are not apparent in the inner (< 1 R e ) stellar population 
gradients studied here, as minor mergers only affect stellar pop¬ 
ulation gradients at > 2 R e ([Foster et al.|2009||La Barbera et ah] 
|20121|Pastorello et al.|2014[|Hirschmann et al.|2014) . However, the 
photometric results are evidence of the active merging activity of 
these galaxies. Furthermore, 4 of the galaxies in our sample have 
close massive companions, most of these companions are FRs and 
are gravitationally bound to their respective BCG. This suggests a 
potential future major merger ( |Jimmy et al.|2013| ). 

Many studies have found similar results on the stellar popu¬ 
lations of BCGs (e.g. |Whiley et al.|2008||Stott et al.|2008| ). |Wen &| 
|Han| ( |20lT) analysed the BCGs colours from different high-redshift 
data sets (CFHT, COSMOS, SWIRE) and found that BCGs are 
consistent with stellar population synthesis models in which the 
galaxy formed at z > 2. However, a large fraction of the sample 
shows bluer colours on the g’ - z’ and B - mbands at z ~ 0.8, 
indicating star formation at those epochs. Furthermore, some 
BCGs show low levels of star formation in the local Universe (e.g. 


Liu et al. 

i2012, Oliva-Altamirano et al. 2014, Fraser-McKelvie 

et al. 2014 

). This is consistent with the hypothesis that BCGs have 


complex accretion scenarios. 


Thanks to the spatial extent of the IFU spectroscopy we were able 
to resolve 3 companion galaxies (1027B, 1048B, 1048C) from 2 of 
the BCGs (1027A, 1048A). We find that the companion galaxies 
similar stellar populations with their respective BCG. However, 
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due to the fact that their effective radius are close to the seeing 
FWHM, their stellar population gradients are unreliable. 

6.4 Connection Between Stellar Populations and Kinematics 

One of the advantages of using IFU spectroscopy for this analysis 
is that we can compare the stellar populations to the kinematics 
of the galaxies in our sample. We have analysed 7 slow and 2 fast 
rotating BCGs and their kinematics appear to be independent of 
their stellar populations. We do not find any correlation between 
the angular momentum of the galaxies and their stellar population 
gradients in this small sample. The SRs show a large scatter in 
their metallicity gradients. The 2 FRs have similar metallicity 
gradients, but within the range of the SRs. 

|Moody et al.| < |2014|) and |Naab et "ah| ( |2014| > showed that if 
the gas fraction in a major merger is less than 10 per cent (dissi¬ 
pationless), the galaxy tends to maintain the slow or fast rotation 
of their progenitors. The BCGs in our sample are likely to have 
preserved the slow or fast rotation of their progenitors, given the 
lack of gas observed at present epochs and the evidence of recent 
dissipationless major mergers. 


7 CONCLUSIONS 

We present the first analysis of the stellar populations of 9 BCGs 
using IFU observations. We compare their stellar populations to 
those of the IFU observations of field early-type galaxies from 
the SAURON and AILAS 3D sample. We draw the following 
conclusions: 

(1) The BCGs have a wide range of stellar ages, high metal- 
licities, and have shallow metallicity gradients. This implies 
diverse evolutionary paths (passive and active accretion). The 
BCGs’ central stellar populations and gradients are consistent with 
those of early-type galaxies of similar mass, with the exception 
that the BCGs have a wide range of central ages. 

(2) Three of the BCGs have similar mass close companions 
galaxies within 18 kpc. From those 3 BCGs we were able to 
resolve the companions of 2 of them. The companion galaxies 
have central stellar populations consistent with their respective 
BCG. 

(3) We do not observe a relationship between the stellar pop¬ 
ulations of BCGs and their stellar kinematics (slow and fast 
rotators). 

IFU analysis has allowed us to determine the angular mo¬ 
mentum of BCGs and to study their stellar population gradients 
without the orientation bias innate to long-slit spectroscopy. It 
has also allowed us to include their companion galaxies in the 
analysis. This study hints of intriguing differences between BCGs 
and similar mass early-type galaxies but requires much larger 
samples to confirm. The SAMI galaxy survey ([Fogarty et al.| 
|2014[ |Allen & the SAMI Galaxy Survey Team||2014| i~ currently 

underway will allow us to achieve this goal. New ultra wide-field 
IFU spectrographs such as MUSE on the Very Large Telescope 
and large date sets as the MASSIVE survey JMa et al.|2014) will 
also allow us to determine whether the picture is different at radii 
beyond 1 R e . 
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APPENDIX A: REPRESENTATIVE LICK INDICES 
EQUIVALENT WIDTH 

In Fig[8]we show the indices H/3, Fe5010, and Mgb, proxies of age 
and metallicity respectively, as a function of velocity dispersion. 
The BCGs show higher H/3 and Fe5010, and similar Mgb values 
than the early-type galaxies at fixed velocity dispersion. 


APPENDIX B: STELLAR POPULATION PROFILES OF 
BRIGHTEST CLUSTER GALAXIES AND THEIR 
COMPANIONS 

In this Section we show the stellar population profiles for each of 
the SR BCGs (Fig[9j, FR BCGs (Fig|TUJ, and FR companion galax¬ 
ies (FigfTT). 
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log G e 


Figure 8. Lick indices as a function of velocity dispersion. The squares represent the BCGs and the blue crosses represent the ATLAS 3£) galaxies. The indices 
were measured at a Lick/IDS resolution (> 8.4 A). The BCGs show higher H/3 and Fe5010, and similar Mgb values than the early-type galaxies at fixed 
velocity dispersion 
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Figure 9. Metallicity (upper panel) and age (lower panel) profiles of the SR BCGs. The name of the galaxy can be found in the upper-left corner. In the 
lower-right corner we show the gradient values. The shaded area represents the seeing FWHM (0.9"). Each green square represents the metallicity and age 
value of each annulus in the galaxy. The dotted line represents the best fit to the profile. The dashed line indicates the central region (an aperture of 0.2 R e ). . 
The SR BCGs have a large scatter in their observed metallicity and age gradients. 
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Figure 10. Metallicity (upper panel) and age (lower panel) profiles of the FR BCGs. The name of the galaxy can be found in the upper-left corner. In the 
lower-right corner we show the gradient values. The dotted line represents the best fit to the profile. The dashed line indicates the central region (an aperture 
of 0.2 R e ). The shaded area represents the seeing FWHM (0.9"). Each red-open square represents the metallicity and age value of each annulus in the galaxy. 
Both of the FR BCGs have negative metallicity gradients and positive age gradients. 
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Figure 11. Metallicity and age profiles of the companion galaxies (all FRs). Each double panel shows the metallicity (upper panel) and age (lower panel) 
profiles of each companion. The name of the galaxy can be found in the upper-left corner. The dotted line represents the best fit to the profile. The dashed line 
delimitates the central region (an aperture of 0.2 R e ). The shaded area represents the seeing FWHM (0.9 // )- Each open star represents the metallicity and age 
value of each annulus in the galaxy. All 3 of them have a flat metallicity gradient. However, in 1048B and 1048C, this could be a result of the seeing FWHM 
being equivalent to the R e in these galaxies. 
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